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The linear instability that induces a relativistic electron beam passing through a return plasma 
current to filament transversely is often related to some filamentation mode with wave vector normal 
to the beam or confused with Weibel modes. We show that these modes may not be relevant in this 
matter and identify the most unstable mode on the two-stream/filamentation branch as the main 
trigger for filamentation. This sets both the characteristic transverse and longitudinal filamentation 
scales in the non-resistive initial stage. 
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Inertial confinement fusion schemes commonly involve 
in their final stage the interaction between some highly 
energetic particle beams and a dense plasma target. 
This is in particular valid for the Fast Ignition Scenario 
1] (FIS) where some laser-produced relativistic electron 
beam would eventually propagate into the dense plasma 
where it would be stopped. This process would lead to 
strong local heating and the ignition of a fusion burn 
wave. In this respect, microscopic turbulence in beam- 
plasma systems is one of the main potentially deleterious 
effects for inertial fusion schemes since it may prevent the 
conditions for burn to be met by broadening the phase 
area where particles deposit their energy. Within the FIS 
framework, a strong research effort has thus been put re- 
cently on the interaction of a relativistic electron beam 
with a plasma with a focus on beam filamentation insta- 
bility, that is microscopic in the transverse direction (see 
e.g. [SUSHI!)- The experimental evidence of filamen- 
tation of very high current laser-produced electron beams 
was recently reported for conditions relevant to the FIS 
■ More generally, filamentation is a potential instabil- 
ity in beam-plasma systems in frameworks ranging from 
accelerators physics to solar flares. 

In the linear stage, filamentation is generally studied 
under some simplifying ab-initio transverse approxima- 
tion of the dielectric tensor, so that filamentation insta- 
bility is attributed to the exponential growth of unstable 
electromagnetic purely transverse modes (k • E = 0) with 
wave vector k normal to the beam 0, H 0> El EH E3 • 
It is also common to refer to this instability as Weibel 
instability S S Eli though the original mode studied 
by Weibel [ljj would require some plasma temperature 
anisotropy to be driven. Figure Q sketches the original 
definitions of various modes under the original Weibel's 
scenario where k is parallel to the beam, along the low 
temperature axis. As long as the beam is not relativis- 
tic, the largest instability it undergoes is the two-stream 
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FIG. 1: Weibel, two-stream and filamentation modes. 



one, where the second "stream" is the return current it 
generates in the plasma. But in the relativistic regime, 
the "filamentation" growth rate eventually exceeds the 
two-stream one and is supposed to induce beam filamen- 
tation. 

In reality, the beam suffers much more instabilities at 
the same time. Indeed, filamentation, Weibel or two- 
stream instabilities pertain to various orientations of the 
wave vector and various kinds of waves (transverse or 
longitudinal) , but in the real world the beam-plasma sys- 
tem triggers every possible modes allowed by Maxwell 
equations with a wide range of wave vector orientation. 
Among all the triggered modes, the unstable ones shall 
start growing exponentially while the most unstable one 
shall mostly shape the beam. When it comes to know 
how the beam is eventually affected when entering the 
plasma, one needs therefore to answer the two questions: 
1) which is the most unstable mode all over the k space 
for the system investigated ? and 2) how shall this mode 
shape the beam ? Following the guideline built by these 
two questions, we assert that the so-called "filamenta- 
tion" instability is not the fastest growing instability, 
even in the relativistic regime, so that it is not the an- 
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swer to the first question. As for the second question, 
we shall see that this instability would not produce the 
observed effects anyway, even if it were the stronger one. 
We shall conclude proposing a new "candidate" for beam 
filamentation and comparing ourpredictions with the ex- 
perimental results presented in Q. For clarity, we shall 
keep labelling the most unstable transverse mode with 
wave vector normal to the beam as the "filamentation" 
mode, though our point is precisely that it does not fila- 
ment. 

Let us consider a beam of electrons (having mass m e 
and charge e) of density n& and relativistic velocity V& 
passing through a return current of plasma electrons of 
density n p , so that the system is unmagnetized. Both 
beam and plasma are infinite and homogenous and ions 
are supposed to form a fixed neutralizing background. 
Let us define the ratio a — n b /n p and introduce the 
plasma frequency uj p = (47rn p e 2 /TO e ) 1 / 2 . Here the beam 
will be assumed to be cold in the longitudinal direc- 
tion, which is correct provided the ratio of its longitu- 
dinal thermal velocity V tb \\ over the parallel phase veloc- 
ity w p /k\\ is small compared to a 1 / 3 . The filamentation 
growth rate can then be evaluated in the weak beam den- 
sity limit (a -C 1) through 
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with (3 = V b /c and j b = - /3 2 . Within the same 

weak beam density limit, the two-stream growth rate 
reads 
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Since 6f decreases like 7 h whereas Sts decreases like 
7 ; T 1 , the filamentation growth rate eventually exceeds the 
two-stream one when the beam is relativistic. Compar- 
ing filamentation growth rate with the Weibel one (trans- 
verse waves with wave vector along the beam, as in |l3|h 
one finds filamentation to be also dominant so that it 
eventually appears to be the largest instability [2l| . 

However, this conclusion needs to be modified when ac- 
counting for every other unstable modes with wave vector 
neither normal nor parallel to the beam. Investigating 
these modes demands a fully electromagnetic formalism 
which is the only way to capture longitudinal modes (two- 
stream) as well as transverse modes (Weibel and filamen- 
tation). Indeed, such a procedure shows that two-stream 
and filamentation modes pertain to the same branch of 
the dispersion equation so that it is possible to switch 
continuously from the former to the later by increasing 
the angle 9k between the beam and the wave vector from 
to 7r/2. Consequently, the angle ipk between the wave 
vector and the electric field of the mode (ipk = (k, E)) 
needs to go continuously from to n/2 to bridge be- 



tween longitudinal two-stream modes and transverse fil- 
amentation modes. In a recent paper [lij ]. we began to 
implement such an electromagnetic formalism using the 
relativistic Vlasov equation to describe the evolution of 
the electronic distribution function of the beam-plasma 
system. Using some simple waterbag distribution func- 
tions for the beam and the plasma, we investigated the 
two-stream/filamentation (TSF) branch and found that 
the growth rate reaches a maximum for an intermediate 
orientation of the wave vector. This maximum scales like 
7 h 1 ^ 3 and reads 
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It is noticeable that this result may be recovered under 
the electrostatic longitudinal approximation [l^ . Such 
an approach cannot however sweep the whole k-plane. 
Equation shows that, even in the relativistic regime, 
filamentation growth rate should not be the larger one. 
On the contrary, this mixed two-stream filamentation 
mode shall be all the more dominant over the usual fil- 
amentation mode that the beam is relativistic because 
of its 7b scaling. This trend amplifies even more when 
accounting for transverse beam temperature, since fila- 
mentation is damped 0,0] while 5m is almost unaffected 
[lfij . Therefore, one can say that the so-called filamenta- 
tion instability may not be the fastest growing one. 

Let us explore this further and move to our second 
point by questioning on what filamentation instability 
would do to the beam, if even it had the largest growth 
rate. Within the linear approximation, one restricts 
to small fluctuations of the electron charge density. If 
Pi and Ei denote the first order perturbations, respec- 
tively, to the electron charge density and to the electric 
field, Poisson equation written in Fourier space brings 
k • Ei(k, lo) — 47rpi(k, lo). It comes directly from this 
equation that a transverse mode with k ■ Ei =0 has 
pi (k, to) = and cannot yield density perturbations 
within the limits of the linear regime. As a consequence, 
the transverse filamentation instability with wave- vector 
normal to the beam cannot yield any charge density fluc- 
tuations from the linear stage. It is important to note 
that, for the same reasons, the original Weibel mode 0] 
cannot linearly induce density filamentation either. 

As far as beam filamentation is concerned, experi- 
ments and simulations show that the electronic density 
varies transversely to the beam, producing the filaments 
@ El ■ Now, if the electronic density varies while 
background ions are (almost) at rest, there is neces- 
sarily a net charge perturbation which precisely cannot 
be accounted for by the mere exponential growth of a 
purely transverse wave. It seems therefore that even if 
it were the fastest growing instability, the so-called fila- 
mentation instability would not produce these filaments. 
It is worth noticing that it could produce current fil- 
aments, for Maxwell's equations allow such a wave to 
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produce such perturbations. But these current filaments 
would have to preserve the neutrality of the system beam- 
plasma, that is to preserve electronic density since ions 
can be considered at rest. 

Let us eventually determine which mode is responsible 
for the observed filamentation. We shall see here that the 
most natural candidate is the most unstable mode found 
along the TSF branch. Being the fastest growing mode, 
it is the one whose growth should "shape" the beam dur- 
ing the linear phase while the other modes create fluc- 
tuations around this basic shape. As for its ability to 
create filaments, it is quasi-longitudinal 0, 0] so that 
its divergence does not vanish. This mode, unlike the so- 
called filamentation mode, satisfies therefore the criteria 
to induce filamentation: It is the fastest growing one, it 
is microscopic in the transverse direction and it is two- 
stream like, that is, quasi-longitudinal. Expressing the 
density perturbation in terms of the wave electric field 
yields pi(k, ui) = kE\ (k, uj) cos ifk/^ir, and one retrieves 
the density perturbation in the real space through 

fcSi(k,u; k )cos^ fc 
p 1 (r,t) = 2_^ j- exp(ik • r - iw^t) (4) 

The sum above runs over every wave vector and every 
proper frequencies Wk- Yet it will obviously be domi- 
nated by the contribution of the fastest growing (hav- 
ing <5k = 3(wk) > 0) self-excited modes. Figure dis- 
plays the growth rates on the TSF branch in the fen) 
plane [22] for some zero or finite plasma thermal velocities 
Vt p — Vtp\\ — Vt p i_ and some zero or finite beam trans- 
verse thermal velocities Vtb± . It is important to note 
that the associated real parts 9?(wk) are in the vicinity of 
the resonance given by u> — kit Vf, = 0. To our knowledge, 
this is the first exact computation of TSF growth rates 
in the whole fc-space including beam and plasma tem- 
peratures effects. These curves clearly show that when 
temperatures are accounted for they act to control the 
instability domain, damping the small wavelengths per- 
turbations along the filamentation direction (fcy = 0) and 
deforming the growth rates surface so that a maximum 
growth rate appears for a finite oblique wave vector kj/ . 
In this respect, Fig. 2(b) shows the drastic influence of 
beam transverse temperature for a cold plasma. Yet, ev- 
ery physical plasma has a finite bulk temperature and, for 
a small enough, this plasma temperature can be shown 
|l6| to control essentially the maximum growth rate lo- 
cation. Its (fcj_,fc||) components are then 

We can then roughly evaluate the density perturbation 
in Eq. (0J by only retaining the k^/ contribution. As 
for the corresponding proper frequency, one has u>k M ~ 
±ujp + iSm where 5m is given by Eq. 1(3} . Summing 
in the k-space the four contributions associated to all 



the possible orientations of kjf in J5J , one finds that the 
density perturbation behaves essentially as 

/9i(r,i) oc exp((5Mi) sin(k M \\z — ui p t) cos(k M ±x). (6) 

Equation © displays spacial modulation of electron den- 
sity in the beam direction (z) as well as in the normal 
direction {%). In the normal direction, we witness the 
"birth" of the beam filamentation in the linear stage, 
with filaments interspace 

L f ~Tr\ a yfVt P /V b , (7) 

where X s = c/u> p is the skin depth. 

There are by now very few relevant experimental re- 
sults available for quantitative comparisons with this re- 
sult. We can consider Fig. 3 of Ref. 0, where plasma 
electronic density is about 10 20 electrons/cm 3 . This 
yields a plasma skin depth of about 53 fim while Fig. 
3 scale indicates the transverse space between filaments 
is somehow smaller. Indeed, the quantity Lf introduced 
above is the skin depth times a Tty/Vtp/Vb factor which is 
smaller than 1 for a non-relativistic plasma since V& ~ c 
here. Taking account of the estimated plasma temper- 
ature (100 eV), we finally find Lf ~ 23 /im which is in 
good agreement with what is observed. Figure displays 
the right hand side of Eq. JJj} for t — l/u> p . Filaments 
are clearly visible, combined with a beam segmentation 
along the beam direction into segments ir\ s Vb/c ~ ir\ s 
long. This parallel segmentation may not be easily dis- 
tinguishable on Fig. 3 of Ref. Q for its characteristic 
length (more than 150 /im) is comparable to the size of 
the entire picture. 

Let us here briefly discuss the applicability of our 
study. As far as plasma (or beam) transverse temper- 
atures are large enough, the above results should apply 
to a finite system where the beam has a finite radial ex- 
tension rb, the condition being that rbk^i± 3> 1 (see 
for filamentation instability in a finite size beam) . As far 
as FIS quantitative applications are concerned, the major 
potential restriction of the present study is the fact that 
the longitudinal beam temperature has been neglected. 
Taking it into account would substantially increase its 
difficulty as it would require a full kinetic treatment and 
may render untractable the already demanding formal 
computations used in Fig. [21 A useful discussion on the 
onset of kinetic effects and breakdown of the cold beam 
hypothesis may be found in Refs. [l^ . Besides, we used 
there waterbag distributions which were simpler to tackle 
than Maxwellian, but this should only affect marginally 
the quantitative results obtained. 

Let us summarize our point as a conclusion. Relativis- 
tic beam filamentation is an observed phenomena. It is 
usually associated with the exponential growth of an un- 
stable mode called "filamentation instability" . It turns 
out that a thorough study of every unstable modes re- 
veals that the "filamentation mode" should not be the 
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FIG. 2: Growth rates in terms of Z = kVb/ui p . (a) Cold beam - cold plasma (see also fT^|L (b) hot beam - cold plasma and (c) 
hot beam - hot plasma. Parameters are a = 0.05 and 7b = 4 for (a,b,c), Vtb± = 14/10 for (b,c) and Vtp = V&/10 for (c). 




FIG. 3: Right hand side of Eq. @ for t — l/u> p . Parameters 
are Vt p = c/30, a = 0.1 and 7 6 = 7 yielding 8m ~ 0.16lj p . 



most unstable. Furthermore, this mode is purely trans- 
verse and therefore unable to produce charge density per- 
turbations. A better candidate to explain beam filamen- 
tation is the most unstable mode all over the k space, 
which turns to be intermediate between filamentation 
and two-stream waves. Not only this mode appears to 
be the fastest growing one, it is also quasi-longitudinal 
so that it can perfectly induce charge density perturba- 
tions. A simple evaluation of its growth shows how it 
creates beam filaments within a few plasma periods and 
agreement with experiment presented in is found to be 
correct. It sets the characteristic transverse and longitu- 
dinal filamentation scales, at least during the linear ini- 
tial stage when resistive (collisional) effects are still neg- 
ligible ,18]. Finally we wish to mention that our study 
emphasizes the importance of quasi-longitudinal modes 
in modelling filamentation which agrees with some con- 
siderations recently put forward by Macchi et al. [13] 
among others. 
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